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We study the conductance of phase-coherent disordered quantum wires focusing on the case 
in which the number of conducting channels is imbalanced between two propagating directions. 
If the number of channels in one direction is by one greater than that in the opposite direction, 
one perfectly conducting channel without backscattering is stabilized regardless of wire length. 
Consequently, the dimensionless conductance does not vanish but converges to unity in the 
long-wire limit, indicating the absence of Anderson localization. To observe the influence of a 
perfectly conducting channel, we numerically obtain the distribution of conductance in both 
cases with and without a perfectly conducting channel. We show that the characteristic form 
of the distribution is notably modified in the presence of a perfectly conducting channel. 
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1. Introduction 

The discovery of a perfectly conducting channel in 
disordered carbon nanotubes^^^ and graphene nanorib- 
bons^'^ provides a counterexample to the conjecture that 
an ordinary quasi-one-dimensional quantum system with 
disorder exhibits Anderson localization. The band struc- 
ture of carbon nanotubes and graphene nanoribbons has 
two energy valleys near the Dirac points, which are well 
separated in momentum space. Ando and co-workers 
have shown that a perfectly conducting channel with- 
out backscattering appears in each valley of a metal- 
lic carbon nanotube when impurity scattering does not 
connect the two valleys.^""' Such a situation can be re- 
alized when only long-ranged impurities are distributed. 
The presence of a perfectly conducting channel is con- 
cluded from the facts that the reflection matrix r for 



each valley has the skew-symmetry (i.e. 



— r) and 



that the number of conducting channels is odd in each 
valley. It has been pointed out that a perfectly conduct- 
ing channel is not intrinsic to carbon nanotubes with 
long-ranged disorder, but universally appears in disor- 
dered quasi-one-dimensional systems which belong to the 
symplectic universality class with an odd number of con- 
ducting channels. ^"^^ Note that the symplectic univer- 
sality class consists of systems having time-reversal sym- 
metry without spin-rotation invariance. The discovery of 
a perfectly conducting channel in the symplectic univer- 
sality class resolves the long-standing puzzle raised by 
Zirnbauer and co-workers. ^^"^^ The presence of a per- 
fectly conducting channel in graphene nanoribbons with 
two zigzag edges has been pointed out by Wakabayashi 
and co-workers.^' ^ As in the case of carbon nanotubes, 
its presence is ensured only when impurities are long- 
ranged. However, the corresponding stabilization mecha- 
nism of a perfectly conducting channel is completely dif- 
ferent from that in the case of carbon nanotubes. In each 
valley of a graphene nanoribbon with two zigzag edges, 
the number of conducting channels in one propagating 



direction is by one greater, or smaller, than that in the 
opposite direction because each valley has an excess one- 
way channel. This imbalance results in the appearance 
of a perfectly conducting channel. The presence of 
a perfectly conducting channel is a universal character- 
istic emerging in disordered quasi-one-dimensional sys- 
tems which belong to the unitary universality class with 
the channel-number imbalance between two propagating 
directions. ^''"^^ The unitary universality class is char- 
acterized by the absence of time-reversal symmetry. 

The behavior of the dimensionless conductance in 
such systems with a perfectly conducting channel 
has been studied by means of the scaling theory 
based on a random-matrix model,^'*'^-'^'^^'^^'^^'^^ the 
super-symmetric field theory,^' the mapping to a 
super-spin chain model, and numerical simula- 
tions. It has been shown that the dimen- 
sionless conductance decays exponentially with increas- 
ing system length L towards unity, indicating the absence 
of Anderson localization. This behavior should be con- 
trasted to that in the ordinary case without a perfectly 
conducting channel, in which the dimensionless conduc- 
tance vanishes in the long-L limit. It has also been shown 
that the dimensionless conductance in the presence of 
a perfectly conducting channel decays much faster than 
that in the ordinary case. To deeply understand such 
peculiar transport properties of the systems with a per- 
fectly conducting channel, it is desirable to study the 
distribution of dimensionless conductance. 

For simplicity, we hereafter restrict our attention to 
disordered quantum wires with unitary symmetry. We 
consider wires of length L with N right-moving channels 
and N -\- m left-moving channels with m = 0, 1. In the 
ordinary case of m = 0, the right-moving channels and 
the left-moving channels have the same set of transmis- 
sion eigenvalues {ri,T2, . . . ,T/v}. Therefore, the dimen- 
sionless conductance Gl and Gr for the left-moving and 
right-moving channels, respectively, are equivalent and 
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are given by Gr = Gl = f) with g = X]o=i In the case 
of m = 1 , one left-moving channel becomes perfectly con- 
ducting and Gl and Gr differ from each other. If the set 
of transmission eigenvalues for the right-moving channels 
is {Ti, T2, . . . , Tjv}, that for the left-moving channels is 
expressed as {Ti, • • • , TV, 1}, where we have identified 
the A^-|-lth channel as the perfectly conducting one.^"^ We 
obtain Gj, = g + 1 and Gr = g with g = J2a=i -^a- "^^^ 
above argument indicates that in both the cases of m = 
and of m = 1, we can completely describe the statistical 
property of Gl and Gr if the distribution of g is given as 
a function of the wire length L. Let p{g) be the distribu- 
tion of g. As the averaged conductance (g) monotonically 
decreases with increasing L, we employ it as the scaling 
parameter instead of L. In the absence of a perfectly con- 
ducting channel (i.e., m = 0), it has been shown that the 
characteristic form of p{g) crucially depends on {g) ; the 
dimensionless conductance g obeys a Gaussian distribu- 
tion in the short-wire regime of (g) ^ 1,^^ while a log- 
normal distribution is expected in the long-wire regime 
of (g) <C 1.^^ It has also been shown that an interesting 
behavior appears in the crossover regime of (g) ~ 1.29,30 
It is interesting how the behavior of p{g) is affected by 
the presence of a perfectly conducting channel. 

In this paper, we theoretically study the behavior of 
p(g) in disordered quantum wires with unitary symme- 
try focusing on the crossover regime of (g) ~ 1. We ob- 
tain p(g) in both the cases of to = and of m = 1 at 
(g) = 1.0, 0.5 and 0.35. To obtain p{g), we first adopt 
a classical Monte Carlo approach based on the existing 
scaling theory. Employing an approximate analytic ex- 
pression of the probability distribution function of trans- 
mission eigenvalues, we calculate p{g) by using a classical 
Monte Carlo method. For comparison, we also calculate 
p{g) by adopting a numerical simulation approach with 
a tight-binding model. For the case of to = 0, we employ 
the tight-binding model for a graphene nanoribbon with 
zigzag and bearded edges, while the tight-binding model 
for a graphene nanoribbon with two zigzag edges is em- 
ployed for the case of to, = 1 . We confirm that the Monte 
Carlo approach and the numerical simulation approach 
provide qualitatively identical results. We show that at 
(g) = 1.0 (i.e., near the diffusive regime), the distribu- 
tion p{g) in the case of m = 1 is similar to that in the 
ordinary case of m = 0. However, a notable difference 
appears for the smaller values of (g). It is shown that in 
the case of to = 0, the distribution p{g) at both (g) = 0.5 
and 0.35 are cut off when g exceeds unity, resulting in the 
appearance of a kink near g = 1,^*^' while no such struc- 
ture appears in the case of to = 1. This indicates that 
the absence of a kink in p{g) is a notable characteristic 
of disordered quantum wires with a perfectly conducting 
channel. We explain this behavior from the viewpoint 
of eigenvalue repulsion due to the perfectly conducting 
transmission eigenvalue. 

In the next section, we obtain the conductance distri- 
bution p{g) by using a classical Monte Carlo approach 
based on the scaling theory. In §3, we obtain p{g) by us- 
ing a numerical simulation approach with a tight-binding 
model. Section 4 is devoted to summary. 



2. Monte Carlo Approach 

We present the scaling theory based on a random- 
matrix model to describe electron transport properties 
in systems with the channel-number imbalance. Let us 
consider disordered wires with N right-moving channels 
and N + m left-moving channels with to = 0, 1. As dis- 
cussed in the previous section, the dimensionless conduc- 
tance Gl and Gr are expressed in terms of g = 
as Gl = m + g and Gr = g. It should be noted that the 
mean free paths I and I' for the left-moving and right- 
moving channels, respectively, are not equal if to ^ 0. 
Indeed, they satisfy V = {N/{N + m))l. The statistical 
behavior of g is described by the probability distribu- 
tion of {Ti,T2,...,Tw}. We define K = {I - Ta)/Ta 
and introduce the probability distribution P{{\a);s) of 
{Ai, A2, . . . , Ajv}, where s is the normalized system length 
defined by s = L/l. The Fokker-Planck equation for 
P({Aa}; s), which is called the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) equation,^^'^^ is expressed as^^ 

dP{{Xa};s) 
ds 

with 

N N-l N 

j=n^"x n n \>^a-x,?- (2) 

c=l 6=1 0=6-1-1 

The factor -^c^ ™ represents the repulsion arising 

from the perfectly conducting eigenvalue. This eigenvalue 
repulsion suppresses non-perfectly conducting eigenval- 
ues {Ti,T2, . . . jTjv}. It should be emphasized that the 
influence of a perfectly conducting channel is described 
by only this factor. 

We numerically calculate p{g) by using a classical 
Monte Carlo approach based on an approximate proba- 
bility distribution of transmission eigenvalues. The prob- 
ability distribution can be obtained from the exact solu- 
tion of the DMPK equation.^^'^** We define Xa by the 
relation of Aa = sinh^ Xa- In terms of a set of variables 
{Xa\, the probability distribution is given by^^ 

P{{xay, s) = const, e-^^i'^-}) (3) 

with 

H{{xa}) = (7^^ " ("^ + |a;asinh2a;a|^ 

(ln|sinh Xa — smh a;6| -|- In |a;^ — a;;, | 1 . 

a,6=l(a>5) 

(4) 

Although eq. (3) is derived under the assumption of 
s/4:N <C 1, we can expect that it is qualitatively reli- 
able even when s/4:N ~ 1?^ Note that g is expressed as 

g = J2a=i l/cosh^Xa. 

We can interpret H{{xa}) as the Hamiltonian function 
of N classical particles in one dimension. This analogy 
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allows us to adapt a Monte Carlo approach to numeri- 
cal calculations of p{g). Using a simple Metropolis algo- 
rithm,^'^' we obtain p{g) for A'' = 5 in both the cases of 
TO = and of TO = 1 at {g) = 1.0, 0.5 and 0.35. By nu- 
merically calculating {g) as a function of the normalized 
system length s, we find that (g) = 1.0 at s/N = 0.77 
(0.45), (g) = 0.5 at s/N = 1.59 (0.75) and (g) = 0.35 at 
s/N = 2.15 (0.91) for TO = (to = 1). As pointed out in 
ref 23, the decrease of (g) as a function of s becomes fast 
with increasing to due to the eigenvalue repulsion arising 
from the perfectly conducting eigenvalue T = 1. There- 
fore, the value of s/N yielding a given (g) for to = 1 is 
smaller than that for to = 0. 



Pis) (a) 

2- 
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Fig. 1. The conductance distributions at (a) (g) = 1.0, (b) ((/) = 
0.5 and (c) (g) = 0.35 for m = (solid lines) and m = 1 (dotted 
lines). 



The resulting conductance distributions are shown in 
Fig. 1. Each distribution is obtained by using 2 x 10^ 
Monte Carlo steps. Except for the case of {g) = 1.0, a 
notable difference appears between the result for to = 
and that for to = 1. In the ordinary case of to = 0, 
the distributions at both {g) = 0.5 and 0.35 are cut off 
when g exceeds unity, resulting in the appearance of a 
kink near g = ]^ 29,30 y^j^ observe that no such structure 
appears when to = 1. The behavior similar to this has 
been observed in the numerical result of refs. 16 and 22. 

We explain this behavior as follows. In the crossover 
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regime, g is dominated by the largest transmission eigen- 
value Ti and the second-largest one T2. That is, g « 
Ti +T2. In the case of to = 0, the largest eigenvalue Ti 
can become nearly equal to unity. However, the second- 
largest eigenvalue T2 is suppressed by the eigenvalue re- 
pulsion from Ti and thus T2 <^ Ti < 1. Therefore, g 
rarely becomes larger than unity, resulting in the ap- 
pearance of the kink near g = 1. In contrast, Ti in the 
case of TO = 1 cannot become nearly equal to unity due 
to the eigenvalue repulsion from the perfectly conduct- 
ing eigenvalue. Therefore, p{g) near g' = 1 is suppressed 
and the kink disappears. We also observe that in the case 
of TO = 1, the distribution is suppressed near g = 0, as 
well as g = 1, and the width of p{g) becomes narrow. 
The suppression near g = simply reflects the fact that 
s/N yielding a given (9) is smaller in the case of to = 1 
than in the case of to = 0. From eq. (4), we understand 
that smaller s/N results in smaller {xa} and therefore 
g tends to become large. Thus, p{g) near g = is more 
suppressed in the case of to = 1 than in the case of to = 0. 

We finally note that in the case of to = without a 
perfectly conducting channel, p{g) at {g) = 0.5 is similar 
to the critical conductance distribution for the quantum 
Hall transition^^ and p{g) at (g) = 0.35 is similar to 
the critical conductance distribution for the Anderson 
transition in three dimensions.^'' 

3. Numerical Simulation Approach 

In this section, we calculate conductance distributions 
by adopting a numerical simulation approach. As model 
systems, we consider a graphene nanoribbon with zigzag 
and bearded edges (zigzag-bearded nanoribbon) and a 
graphene nanoribbon with two zigzag edges (zigzag- 
zigzag nanoribbon) (see Fig. 2). In the following, we fix 




(b) 



Fig. 2. The structures of (a) zigzag-bearded nanoribbon and (b) 
zigzag-zigzag nanoribbon. 



the number of lattice sites in the transverse direction 
to be M = 30. We describe electronic states of these 
nanoribbons by using a tight-binding model on the hon- 
eycomb lattice. For the case of to = 0, we employ the 
tight-binding model for the zigzag-bearded nanoribbon, 
while that for the zigzag-zigzag nanoribbon is employed 
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for the case of m = 1. The tight-binding model is given 
by 

H=-J2j,j\i){j\+^V,\i){i\, (5) 

where \i) and Vi represent the localized electron state and 
the impurity potential, respectively, on site i, and jij is 
the transfer integral between sites i and j. We assume 
that 7i J = t ifi and j are nearest neighbors and jij = 
otherwise. The band structures obtained from this tight- 
binding model without disorder are displayed in Fig. 3. 
We observe that the zigzag-bearded nanoribbon has a 
completely flat band,^'^^ while a partially flat band ap- 
pears in the zigzag-zigzag nanoribbon, and that both 
of them commonly possess two energy valleys. The right 
and left valleys are referred to as and K- valleys, 
respectively. We assume that impurity potential arising 
from each scatterer has a spatial range greater than the 
lattice constant a. In this case, inter-valley scattering is 
not induced by such long-ranged impurity potential, and 
the two energy valleys are disconnected. Thus, these en- 
ergy valleys independently contribute electron transport. 
Hereafter, we restrict our consideration to the case in 
which the Fermi energy E is placed within < E/t < 1. 
Note that in the zigzag-zigzag nanoribbon, the partially 
flat band provides an excess left-moving (right-moving) 
channel in the (A-) valley. Therefore, the number 
A^L of left-moving channels is by one greater (smaller) 
than the number TYr of right-moving channels in the 
(i^_) valley. That is, TVl = A^r + 1 in the K+ valley and 




Fig. 3. The band structures of (a) zigzag-bearded nanoribbon 
and (b) zigzag-zigzag nanoribbon with M = 30. 
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+ 1 = .^R in the K- valley. This results in the appear- 
ance of one perfectly conducting channel in both the left- 
moving channels in the valley and the right-moving 
channels in the K- valley.^' ^ In contrast, A^l and Nn 
are equal to each other in both the valleys of the zigzag- 
bearded nanoribbon and no perfectly conducting channel 
appears. Let and be the dimcnsionless conduc- 
tance for the left-moving channels and that for the right- 
moving channels in the valley K±. For the case of the 
zigzag-bearded nanoribbon with A^l = -^R = A^ in both 
the valleys, we observe that G^ ~ ~ G^ — = g 
with g = J2a=i ^a- In contrast, for the case of the zigzag- 
zigzag nanoribbon with A^l — 1 = A^r = N in the 
valley and A''l = — 1 = A'' in the K- valley (see Table 
I) , we observe that G^ = G£ = g and G^ = Gj^ = 1 + g 
with g = J2a=i ^o- Note that G^ and G^ contains 
the contribution from one perfectly conducting channel. 
The total conductance of these systems is expressed as 
Gtotai = m + 2g with m = for the zigzag-bearded 
nanoribbon and m = 1 for the zigzag-zigzag nanoribbon. 



Table I. The number of left- moving channels, ATl , and the num- 
ber of right-moving channels, Wr, where m = for the zigzag- 
bearded nanoribbon and m = 1 for the zigzag-zigzag nanoribbon. 
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K+ 


N + m 


N 




N 


N + m 



Let us consider nanoribbons which is infinitely long 
in the longitudinal direction. To introduce disorder, we 
randomly distribute scatterers in the finite region of L 
lattice sites. We identify L with the length of our system. 
We assume that the impurity potential on site i arising 
from a scatterer at site j is given by 

Fi(i)=«;,-exp(-|n-r,f/rf2), (6) 

where rj is the position vector of site i, and Wj and d 

represent the amplitude and the spatial range of this po- 
tential, respectively. Thus, the impurity potential on site 
i is given by 

Vi^Y.^^{j\ (7) 
3 

where j is summed over the disordered region of length 
L. The amplitude Wj is distributed uniformly within the 
range of —W/2 < Wj < W/2. Let p be the probability 
that each site is occupied by a scatterer. The strength of 
disorder is controlled by W and p. Adopting the model 
described above, we calculate the dimcnsionless conduc- 
tance Gtotai by using the Landauer formula 

Gtotai = tr{ttt}, (8) 

where t is the transmission matrix through the disor- 
dered region. The dimensions of t are equal to {2N+m) x 
{2N+m). We consider the behavior ofg = {Gtoted—Tn)/2 
in the case of A'' = 5. With M = 30, this is realized when 
0.486 < E/t < 0.571 for the zigzag-bearded nanoribbon 
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(m = 0) and 0.529 < E/t < 0.611 for the zigzag-zigzag 
nanoribbon (m — 1). The following parameters are em- 
ployed: W/t = 0.1, p = 0.1 and d/a =1.5 with a being 
the lattice constant (see Fig 2(b)). We can numerically 
obtain t for a given impurity configuration by using a 
recursive Green's function technique. 
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Fig. 4. The average of g for N = 5 with (a) rn = and (b) m = 1 
as a function of the length L, where E/t = 0.49, 0.50, 0.51 and 
0.52 from bottom to top in the case of m = and E/t = 0.54, 
0.55, 0.56 and 0.57 from bottom to top in the case of m = 1. 
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Fig. 5. The average of In 3 for N = 5 with (a) m = and (b) 
m = 1 as a function of the length L, where E/t = 0.49, 0.50, 
0.51 and 0.52 from bottom to top in the case of m = and 
E/t = 0.54, 0.55, 0.56 and 0.57 from bottom to top in the case 
of m = 1. 



Before studying the conductance distribution, wc con- 
sider the behavior of the ensemble averages of g and In g. 



We numerically calculate (g) and (In g) as a function of L 
at E/t = 0.49, 0.50, 0.51 and 0.52 for the case of m = 
and at E/t = 0.54, 0.55, 0.56 and 0.57 for the case of 
m = 1. Figures 4 and 5 show the L-dependence of {g) 
and (In 5), respectively. The ensemble average for each 
data point is performed over 5000 samples with different 
impurity configuration. We observe from Fig. 4 that the 
decay of (g) with increasing L becomes slower with in- 
creasing E/t m both the cases of rn = and of m = 1. 
This implies that the mean free path / becomes longer 
with increasing E/t. Although the decay of {g) seems to 
be slightly faster in the case of rn = 1 than in the or- 
dinary case of m = 0, we cannot state this definitively. 
However, Fig. 5 clearly shows that (In g) decreases faster 
in the case of m = 1 than in the case of m = 0. This is 
consistent with the scaling theory^'^ and is explained as 
follows. In the case of m = 1, one perfectly conducting 
channel is present and the corresponding transmission 
eigenvalue (i.e., T = 1) suppresses the other transmission 
eigenvalues {Ti, T2, . . . , T^q} due to eigenvalue repulsion. 
This results in the suppression of g. Indeed, the scaling 
theory predicts that the conductance decay length ^ de- 
fined by exp[(ln£()] = exp[— 2L/^] is given by^^ 

i = ^ (9) 
^ m+1' ^' 

indicating that ^ in the case of m = 1 is twice shorter 

than that in the case of m = if the mean free path is 
common to the two cases. We estimate ^ from the result 
of (Ing"). The estimated values of ^ are listed in Table II. 
Because the mean free path should depend on E/t and m, 
a detailed comparison with eq. (9) is not straightforward. 
However, we observe the obvious tendency that ^ in the 
case of m = 1 is shorter than that in the case of m = 
0. This is consistent with the prediction of the scaling 
theory. 



Table II. The conductance decay length g in the cases of m = 
and of m = 1 for several values of E/t. 



m 


= 


m 


= 1 


E/t 




E/t 




0.49 


11500 


0.54 


6400 


0.50 


13500 


0.55 


8400 


0.51 


15100 


0.56 


9500 


0.52 


15900 


0.57 


10600 



Now, we consider the conductance distribution p{g). 
We numerically obtain p{g) at (g) = 1.0, 0.5 and 0.35 in 
both the cases of m = and of m = 1. Wc set E/t = 0.49 
for m = and E/t = 0.54 for m = 1. From the numerical 
result displayed in Fig. 4, we find that (g) = 1.0 at L = 
2511 (2583), {g) = 0.5 at L = 5352 (4623) and (g) = 0.35 
at L = 7312 (5801) for m = (m = 1). The resulting 
conductance distributions for the case of m = and those 
for the case of rn = 1 are displayed in Figs. 6 and 7, 
respectively, where each distribution is obtained by using 
5 X 10^ samples. For comparison, we also display the 
conductance distributions obtained by the Monte Carlo 
approach. We observe that the result of the numerical 
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Fig. 6. The conductance distributions for the case of m = at (a) 
(g) = 1.0, (b) (g) = 0.5 and (c) (g) = 0.35. Open circles represent 
the numerical simulation result, while solid lines represent the 
Monte Carlo result. 



Fig. 7. The conductance distributions for the case of rra = 1 at (a) 
{g) = 1.0, (b) (g) = 0.5 and (c) (g) = 0.35. Open circles represent 
the numerical simulation result, while solid lines represent the 
Monte Carlo result. 



simulation approach is qualitatively identical to that of 
the Monte Carlo approach for both the cases of to = 
and of m = 1. Particularly, we again observe a kink near 
g = 1 in the case of m = at (g) = 0.5 and 0.35, while 
no such structure appears in the case of to = 1. 

We finally note that a relatively large deviation be- 
tween the two results appears in the case of m = 1 for 
(g) = 0.5 and 0.35, while no such deviation is observed 
in the case of to = 0. This reflects the fact that the 
probability distribution P({xa}; s) given in eq. (3) over- 
estimates the influence of a perfectly conducting channel 
in the long- wire regime of s/N > 1^^ although this is a 
very good approximation irrespective to s/N in the case 
of m = without a perfectly conducting channel.^" 

4. Summary 

The dimensionless conductance g in disordered quan- 
tum wires with unitary symmetry is studied for both 
cases with and without a perfectly conducting channel. 
To observe the influence of a perfectly conducting chan- 
nel on the behavior of g, we have calculated the con- 
ductance distribution p{g) for these two cases in the 
crossover regime where (g) is slightly smaller than unity. 
We have adopted two approaches, a classical Monte Carlo 
approach based on the existing scaling theory and a nu- 



merical simulation approach with a tight-binding model. 
In the latter approach, we have employed the tight- 
binding model for a graphene nanoribbon with two zigzag 
edges (zigzag and bearded edges) for the case with (with- 
out) a perfectly conducting channel. We have confirmed 
that the two approaches provide qualitatively identical 
resiilts. It is shown that in the absence of a perfectly con- 
ducting channel, the distribution p{g) are cut off when g 
exceeds unity, resulting in the appearance of a kink near 
g = I, while no such structure appears in the presence of 
a perfectly conducting channel. This indicates that the 
absence of a kink in p{g) is a notable characteristic of 
disordered quantum wires with a perfectly conducting 
channel. 
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